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ABSTRACT 

We consider overlap-Dirac fermions at non-zero bare coupling and 
for a small hopping parameter, or, equivalently, large \M\ with M the 
domain- wall height. We prove the existence of a phase at large positive M 
where the abelian axial group {7a (1) is a symmetry, and the corresponding 
pseudo-scalar is an exact Goldstone boson. We also provide a conjecture 
for the phase diagram of asymptotically free gauge theories with overlap- 
Dirac fermions. In particular, we suggest that, for large gauge coupling, 
the massive-fermion phase at negative M possibly extends to all M < 4. 



* Permanent address 



1. Introduction. In recent years it has become clear that the chiral limit can be 
separated from the continuum limit in QCD-like lattice gauge theories [1-15]. Domain- 
wall fermions [1-3] provide a simple way of approaching this limit, and have become 
popular in numerical simulations (see ref. [|J for a review). On the theoretical side, 
a key role in the formulation of the chiral limit is played by the algebraic Ginsparg- 
Wilson relation H (see ref. || for a review). An extensively studied solution of the 
Ginsparg- Wilson relation is the overlap-Dirac operator |7], |J, which may be regarded 
as a certain limiting case of domain- wall fermions |l(| . 

A common feature of all numerical simulations which attempt to approach the 
chiral limit is that they show a pattern of increasing chiral symmetry violations when 
the bare coupling is increased. Keeping such symmetry violations under control in the 
most economic way is important, and analytic results valid beyond the weak-coupling 
limit can provide valuable clues as to what is the best way of achieving this goal. 

In this paper we consider a euclidean SU (N c ) lattice gauge theory with iVj copies 
of overlap-Dirac fermions and with a finite bare coupling g. For a small hopping 
parameter k, the theory may be in a massive or in a massless phase (depending 
on whether A — 1 or A — — 1 in eq. @ below). Our main result is that, in the 
massless phase, the theory may be reformulated such that the abelian axial group 
£7a(1) becomes a manifest symmetry of both the action and the measure, without 
spoiling the locality of the theory. As a result, there are Nf lattice pions in that 
massless phase. The term lattice pion here denotes a pseudo-scalar state which is an 
exact Goldstone boson for finite g. 

To avoid additional complications which are besides our main point, we mainly 
consider lattice theories where the limit g — > defines a confining continuum theory. 
We then expect chiral symmetry to be spontaneously broken in the continuum limit as 
well |1| . (Basically this means that the number of flavors is small enough compared 
to N c ; see below for a comment on the case N c = 3.) If the continuum theory has Nf 
flavors, the number of pions is Nj — 1. The small- K phase referred to above thus has 
one extra (lattice) pion compared to QCD with Nf = iVj. This indicates some sort of 
species doubling, and we will argue that, if the continuum limit is taken inside that 
phase, the number of flavors is actually Nf = 1QN, as for naive fermions. 



Our work was motivated by two recent strong-coupling calculations \ Tf[ fL8j. The 
results of those papers will be compared with ours in the last section. 

This paper is organized as follows. In Sect. 2 we review some relevant properties 
of the overlap-Dirac operator, discussing in particular the physical significance of 



Luscher's axial transformations [D], In Sect. 3 we prove the main result, and in 
Sect. 4 we comment on the role of the index of the overlap-Dirac operator. In Sect. 5 
we discuss the continuum limit in some more detail. In Sect. 6 we turn the existing 
results into a conjecture on the phase diagram of overlap fermions. Finally, in Sect. 7 
we compare our results with those of refs. TB\, and list some issues for future 
research. Locality of the reformulated action is proved in the Appendix. 
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2. The overlap-Dime operator. The Ginsparg- Wilson relation is (we work in units of 
the lattice spacing) M 

D l5 + l5 D = D l5 D. (1) 

The overlap-Dirac operator || which satisfies the Ginsparg- Wilson relation is defined 
by 

D = l + 7s H/\H\, (2) 



(3) 



where X is the usual Wilson-Dirac operator, and 



{Cyi)x,y — n {^x+fi,yU x ,^ $x-[i,yUy „) , (4) 



Bx,y — 9 ^2(^x+fi,yU Xjl _ l + Sx-^yUy ^ . (5) 

Since D is unchanged if H is multiplied by an arbitrary positive number, we use this 
freedom to set A = ±1 in eq. The hopping parameter k is positive by convention. 
(A sign flip of k can be undone by the transformation ip x — ► (— i^i+^a+^s+a* ip x an d 
similarly for if) x .) We take the Wilson parameter to be r = 1, but the discussion can 
easily be generalized to other values. 

The action S(ip, %p) = $ D ip is invariant under Liischer's gauge-field dependent 
axial transformation |TT|, |6| 

8^ = T%^, % = l5 (l-D) = -H/\H\, 

Sip = ipT%. (6) 

Here T is a U(Ni) generator acting on the flavor indices. Since (7s) 2 = 1, Liischer's ax- 
ial transformations together with the usual vector transformations generate a Ul(Ni) x 
U r (Ni) = U v (l) x U A {1) x SU L (N t ) x SU R (Ni) symmetry of the lattice action. This 
statement is true for any k. To avoid confusion we recall that if k = 0(1) the fermion 
measure is in general not invariant under Liischer's Ua(1), leading to the expected 
axial anomaly [|ll|] and a massive rj' particle. 



The symmetry (y) does not always have the physical significance of an axial 
symmetry. Let us first recall the quark spectrum described by a single overlap-Dirac 
field. As usual, the quark spectrum is determined by the free theory. We then have 
for Bq and C , in momentum space, 

Co M (p) = * sin Vp. , B Q (p) = J2 cos Pa* • ( 7 ) 

It is easy to check that at the corners p c of the Brillouin zone (p Cfl = or tt, all /x), 
and only there, D (p c ) is equal to either or 2. A given corner of the Brillouin zone 
gives rise to a massless quark field in the continuum limit if and only if D (p c ) = 0. 



3 



In the context of domain-wall or overlap-Dirac fermions, the customary parame- 
trization of the Wilson-Dirac operator in eq. (Q) is 

X = 4-M + Y,1,C,-B. (8) 

Comparing eqs. (|f) and we see that k = 1/|4 — M\, and that A = 1 (A = — 1) 
corresponds to M < 4 (M > 4). When n components of the four- momentum are equal 
to 7i and the rest are zero, one has Bo(p c ) = 4 — 2n. Referring to the parametrization 
of eq. (H), it follows that for M < all corners have Do(p c ) = 2, while for M > 8 
all corners have Do(p c ) = 0. For M < there are no massless quarks, and when 
M is increased above the values 0, 2, 4, 6 and 8, the numbers of massless quarks that 
are added to the spectrum are 1,4,6,4 and 1 respectively. For M > 8 there are 16 
massless quarks. Note that the points M = 0,2, ... ,8, represent discontinuities in 
the spectrum. 

Suppose now that D (p c ) = at some corner of the Brillouin zone. Near that 
corner Luscher's transformation @ reduces to an ordinary axial transformation (cf. 
eq. (HID below) when acting on the corresponding massless-quark state. We recall that 
the physical axial charge is equal or opposite to the lattice axial charge depending on 
whether n is even or odd, respectively [f[9| . 

The other possibility is that D (p c ) = 2 at a corner. In this case, eq. fl6|) reduces 
to Sip = — T75 ip, Sip = ip T75. Because of the minus sign in the ip transformation rule, 
this is no longer an axial transformation. There exists another version of Luscher's 
transformations where the if) and if) rules look more symmetric, given by (Tl| dip 



T7 5 (l — ^D)ip, Sip = ip T(l — |D)7 5 . In this form, the transformation still reduces to 
an ordinary axial transformation when D (p c ) = 0. But for D (p c ) = 2 this becomes 
Sip = Sip = 0. In other words, for Do(p c ) = 2 the transformation does not act at all 
on states with a momentum close to p c . 

One can summarize the situation by saying that, unlike an ordinary axial sym- 
metry, Luscher's symmetry by itself does not imply the existence of massless quarks. 
But if massless quarks exist, it acts on them as an ordinary axial symmetry. In fact, 
these requirements single out the transformations (|5]) almost uniquely. For consider 
taking Sip = T7 5 (l — 0)ip for some O, with O = for p = 0, but possibly non-zero 
for other momenta in order to accommodate the removal of doublers. We also want 
to form the same Lie algebra as in the continuum, so we need (75(1 — O) ) 2 = 1. From 
this it follows immediately that O satisfies the Ginsparg- Wilson relation. 

3. The massless small-n phase. The properties of D for a small hopping parameter 
k rest on the fact that X = ±1 + O(k) (cf. eq. (|3|)). It follows immediately |Tj]] that 
D = 2 — O(k) for A = 1 while D = O(n) for A = — 1 (these statements refer to the 
eigenvalues of D; we recall that \D\ < 2 always, see ref. ||). In the case A = 1, since 
D = 2 — 0(k) regardless of g, all correlation lengths in the fermion sector must be 
finite in lattice units, and moreover tend to zero for k — > 0. Hence A = 1 corresponds 
to a massive-fermion phase. (This is the usual situation in a hopping expansion. It 
can be established using the methods of ref. pO , or, more explicitly, by applying the 
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techniques of the Appendix to obtain bounds on the kernel of D~ l for A = 1 and 
small k.) 

In this section we will show that A = — 1 corresponds to a massless phase with 
Nf lattice pions. The first important observation |T7) is that the leading, O(k), 



term in the expansion of D is proportional to the naive-fermion operator. Also, we 
have already seen that for small k (large M) overlap-Dirac fermions undergo the same, 
maximal, doubling in the continuum limit as naive fermions. (However, overlap-Dirac 
and naive fermions have a different massless spectrum (of gauge invariant states) at 
finite k and g, see below. Note that the k — » limit is singular for normalized 
expectation values: in this limit there are massless fermions, whereas setting k = 
gives rise to no propagation at all. The hopping expansion in the A = — 1 phase is 
thus qualitatively different from that for Wilson fermions, where, as in the A = 1 
phase, it is an expansion around an infinite fermion-mass theory.) 

In order to analyze the situation for a small but finite hopping parameter it is 



convenient to introduce new variables 17 



9 =(2-D)V, ? = (9) 

In terms of the new variables, Luscher's transformation (§) reduces to the ordinary 
axial transformation 

5q = T'j 5 q, Sq = qT'j 5 . (10) 
The partition function is rewritten as 

Z = det(D) = Jv^Vi) exp(-S(V>,^)) 

= J VqVq det(2 - D) exp(-S'(q,q)) , (11) 

where 

S'(q,q)=J2qD(2-D)- 1 q. (12) 

It is easily verified that D(2 — D)^ 1 anticommutes with 75, as it must, in view of 
eq. (0). The action S'(q,q) is thus invariant under ordinary axial transformations, 
as well as under the usual hyper-cubic rotations. 

The transition to the new variables cannot be done for arbitrary k and A. When k 
exceeds some critical value, (2 — D)^ 1 becomes singular, and in the free fermion limit 
there are poles in the action S'(q,q). Moreover, in non-trivial topological sectors 
S'(q, q) is undefined since there exist stable eigenmodes with eigenvalues and 2. 
Finally, when A — 1, poles appear in S'(q,q) already in the hopping expansion, and 
therefore the discussion below is not applicable in the massive phase. 

When A — — 1, the transition to the new variables turns out to be a powerful 
tool, since in this case eq. (IT3) is a local action for small enough k. In more detail, 
for a range k < kq of the hopping parameter D(2 — D)^ 1 is bounded and has an 
exponentially decaying kernel (as is the standard practice in this context, the last 
statement defines the notion of locality used in this paper). Boundedness is obvious 
since D = O(k), while locality of D(2 — D)~ l is established in the Appendix for 
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k < k (we note that the "admissibility" constraint on the gauge field |L3| is not 
necessary for small k). As for the factor det(2 — D) in eq. QTTD, it is not expected 
to change the universality class (an argument similar to that in the appendix shows 
that trlog(2 — D) is local). 

Thus, we find that for k < k the theory can be consistently formulated in terms 
of a local action S'(q,q), where the action and, obviously, the measure are invariant 
under ordinary vector and axial transformations. The fact that we are dealing here 
with the simple axial transformation of eq. (|H]) is important. As mentioned in the 
introduction, we restrict the discussion to those cases where the limit g —>■ defines a 
confining continuum theory. We then have confinement for any g, and the standard 



arguments that confinement implies chiral symmetry breaking apply []IE |. In the 
present context, the formation of a (qq) condensate was confirmed in the strong 
coupling limit in a 1/N C expansion |T7J (see also ref. jpl] ). 

Therefore, the lattice theory has Nf axial generators, and the pseudo-scalar Gold- 
stone bosons must be in one-to-one correspondence with those. We thus conclude that 
there are Nf lattice pions for k < kq and any (finite) g. The (confining) small- k mass- 
less phase may actually be defined as the part of the phase diagram with Nf lattice 
pions. We expect this phase to extend beyond the region n < n and, in fact, beyond 
the region where (2 — D)~ x is bounded, see below. We comment in passing that in the 



case of naive fermions the lattice symmetry is bigger |]22| , and hence the number of 
lattice pions is larger as well. This naive-fermion symmetry requires that the action 
has only odd-neighbor couplings. 

In summary, we have established the existence of a class of overlap-Dirac theories 
where for A = — 1 and small k there is only one phase for all values of g. This phase 
is characterized by confinement and chiral symmetry breaking, and has Nf lattice 
pions. 

4- The index and Ua{^-)- The question of whether Ua{1) is a symmetry or not may 
be approached from another direction. The integrated Ward identity associated with 



a UaW transformation reads 11 



(SO) = (Tr(%) O) , (13) 

where O = 0(ip,tfj,U^) and the axial variation denoted by 5 is defined in eq. (|G]). 
(To make the above statement well defined we may assume that we work in a finite 
volume.) It is easy to see that Ty{%) = for small n for all gauge field configurations, 
regardless of the sign of A in eq. (|j). First, Tr(7 5 ) is (proportional to) the index of 



D [ITTJ. Then, a non-zero index requires the simultaneous existence of eigenvectors 
with eigenvalues zero and two |6|, [| |14[. But, for A = — 1, D = 0(k) and there can 
be no eigenvalue equal to two. For A — 1, D = 2 — 0(k) and there can be no zero 
eigenvalue. 

What we learn is that the small-K global symmetry of both the A = 1 and the 
A = — 1 phases is Ul(Ni) x Ur(Ni). However, as already explained above, there is 
an important difference between the two phases. In the A = — 1 phase the conserved 
axial generators may be taken to be the ordinary axial generators associated with 



6 



eq. (PI) . We have seen that this implies the existence of massless quarks and (if there 
is confinement) spontaneous chiral symmetry breaking. In contrast, in the A = 1 
phase we have only Liischer's symmetry at our disposal, and as explained in Sect. 2, 
this is not incompatible with the absence of massless-fermion states. 

5. Relation to the continuum limit. The previous results provide valuable information 
for the task of mapping out the (M, g) phase diagram of overlap-Dirac fermions. As 
already mentioned, the hopping expansion is an expansion in 1/|4 — M\. The massless 
(A = —1) phase at small k corresponds to M > M c+ (g), while the massive (A = 1) 
phase corresponds to M < M c -(g). 

We now wish to determine the end points of these two critical lines in the con- 
tinuum limit. We will argue that the massless phase at large positive M extends 
down to M = 8 for g — > 0, whereas the (massive-fermion) phase at large negative M 
extends up to M = provided Ni > 2. The region 0(g) < M < 8 — 0(g) at small g 
is filled with a phase supporting Nf — 1 lattice pions. 

Let us begin with the last statement. For < M < 8 there are massless quarks, 
and since Liischer's symmetry (|6]) acts on those as an axial symmetry, spontaneous 
symmetry breaking should occur, and the corresponding Goldstone bosons should 
exist. In this range no massless state in the Ua(1) channel is expected [|3|, [H| because 



the index of D can be non-zero. Since there are Nf — 1 conserved axial generators 
for non-zero g, this must also be the number of lattice pions. A corollary is that the 
massless large-M phase with Nf lattice pions cannot end at any M < 8 for g — > 0. 
Recall that, according to the free- field (or weak-coupling) analysis, the largest value 
where a phase transition takes place in the g — > limit is M = 8. Assuming that 
that analysis exhausts all possible phase transition points for g — > 0, it follows that 
the large-M phase must end at M — 8. Similarly, the massive-fermion phase with 
no lattice pions at negative M (which is actually a pure glue phase) should end at 
M = for g — > 0. The last statement is true except for A; = 1, where both the 
negative-M and intermediate-M regions do not support lattice pions, and therefore 
they may be analytically connected at non-zero g. 

For small g and M > 8 there are massless states at all the corners of the Brillouin 
zone, namely 16 quark fields per each overlap-Dirac fermion. The lattice U^(l) sym- 



metry corresponds to a flavor non-diagonal axial symmetry in the continuum limit p9 
(see Sect. 2). The total number of pions in the continuum limit is (16 A;) 2 — 1, whereas 
for finite g there are only Nf pions. The difference is explained by states which are 
approximate Goldstone bosons for small g, and which become massless only in the 
continuum limit. This bears some resemblance to staggered fermions, where only one 
pion is exactly massless at finite lattice spacing, while the rest become massless only 
in the continuum limit. 

If one relaxes the assumption we have made in the introduction about the fermion 
content, one can find cases where asymptotic freedom is lost and/or there is no chiral 
symmetry breaking for small g. In such cases there has to be (at least) one phase 
transition as a function of g even at large (positive) M. For the physical case A c = 3, 
asymptotic freedom at large M is lost for N = 2 (Nf = 32), and it is not clear if 
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there is (confinement and) chiral symmetry breaking even for iVj = 1 (Nf = 16) in the 
continuum limit. At strong coupling one always has confinement, and we expect our 
conclusion to hold on the large-g side of such a transition. Also, whatever replaces 
confinement and chiral symmetry breaking at small g should still be consistent with 
having 16AT ; 2 massless fermions for g = 0. Note that this issue does not affect the 
phase diagram in the region of interest for Lattice QCD, with < M < 2. 

6. Conjectured phase diagram of overlap-Dirac fermions. In this section we put 
forward a conjecture for the (M, g) phase diagram of overlap-Dirac fermions. We 
have already discussed three out of four boundaries of the phase diagram. We invoke 
a mean- field argument to cover the last, g — oo, boundary, and then consider the 
simplest finite- (M, g) interpolation. 



oo 




Figure 1: Conjectured phase diagram of overlap-Dirac fermions for Ni > 2. The 
number of lattice pions in each phase is indicated. 

The conjectured phase diagram is depicted in Fig. 1 for JVj > 2. It contains 
the two phases with zero and Nf lattice pions whose existence was established for 
large-negative and large-positive M respectively. There is also a single phase at 
intermediate values of M, characterized by the existence of Nf — 1 lattice pions. The 
phase transition lines emanating from the points M = 2,4,6 must be there because 
the quark spectrum in the continuum limit changes at those points. As explained 
earlier, the regions surrounding these points at small g > all support Nf — 1 lattice 
pions. In our view the most plausible scenario is that these regions are analytically 
connected at non-zero g, but this does not have to be the case. The phase transitions 
are discontinuous at g — 0, and therefore they are discontinuous for small non-zero g, 
by continuity. We expect them to be discontinuous for all g. We remind the reader 
that this phase diagram is applicable provided there is a confining continuum theory 
in the limit g — > for all M. 

If D satisfies the Ginsparg- Wilson relation, so does 2 — D. The passage from D 
to 2 — D is effected by X — > —X, and the phase diagram for 2 — D is obtained by 
the replacement M — > 8 — M. (In the D-phase diagram, however, the critical lines 
may not transform into each other under M — > 8 — M, and in that sense Fig. 1 may 
be misleading.) 
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The symmetries of the three phases are the following. (Below G — > H denotes that 
the global symmetry is G and that H is the symmetry of the vacuum.) The phase 
with Nf lattice pions is characterized by Ul(N{) x Ur(Ni) Uy(Ni). The phase 
with Nf - 1 lattice pion is characterized by U v (l) x SU L (N t ) x SU R (N{) -> U V (N{). 
Finally, the phase with no lattice pions has the global symmetry Ul(N{) x Ur(N{) 
but no massless fermions, and no spontaneous symmetry breaking. 

As can be seen from Fig. 1, we believe that the two small-K phases meet at M = 4 
for g = oo. If this is true, the phase with Nf — 1 lattice pion does not extend to 
g = oo. (It could be that the triple point is at g = oo. Note that one is usually 
interested in taking the continuum limit inside the middle phase because only in the 
interval < M < 2 there is one quark per one overlap fermion in the continuum 
limit.) Let us define (U Xfl ) in a gauge-invariant way, for instance as the fourth root of 
the average plaquette. Then we have that (U Xfl ) — > for g — > oo and, if a mean-field 
analysis is reliable, the g — > oo limit should behave like k — > 0. For M > 4 this means 
that we are in the massless small-K phase with A = —1, which has Nf lattice pions, 
and for M < 4 we are in the A = 1 massive phase, with no lattice pions. 

If correct, the phase diagram of Fig. 1 leads to an important observation. For 
< M < 4, as we increase g at fixed M, we eventually move into the massive phase, 
thus loosing all the massless quarks (even though the action is still invariant under 
Liischer's symmetry)! Remembering the connection with domain- wall fermions ||10|| , 
this means that under the same conditions domain-wall fermions will support (no 
massless quarks and) no massless pions even if the extent of the fifth dimension is 
taken to be arbitrarily large. 

The phase diagram of Wilson fermions is known to contain the Aoki phase [^4| . 



where parity and vector-like symmetries are broken by a flavor non-singlet pseudo- 
scalar condensate. The question may arise whether a similar phenomenon can take 
place in the present context as well. The answer is no. As explained in ref. | f25f , 
the Aoki phase is related to those terms in the effective chiral lagrangian (for finite 
lattice spacing) that break axial symmetries explicitly. There are no such terms in 
the present situation since one has a chiral symmetry at finite lattice spacing as 
well. Using the lattice chiral symmetries one can rotate any non-singlet pseudo-scalar 
condensate into the usual singlet scalar one (qq). An Aoki phase may, however, arise if 
an explicit quark mass term is added to the action. This was investigated numerically 
in the context of domain-wall fermions in ref. 



7. Discussion. Finally, we compare our results with previous work. Our results are 
in agreement with those of Ichinose and Nagao who studied the massless phase to 
second order in the hopping-parameter expansion |T7|. The fate of the Ua(1) pseudo- 



scalar state was left open in ref. [17|, and it was conjectured that it may eventually 
pick up a non-zero mass (in a higher order in 1 / N c ) . As we showed rigorously in this 
paper, in fact U^(l) is an exact symmetry for a small hopping parameter, and the 
corresponding pseudo-scalar is an exact Goldstone boson for any value of g. This 
result is consistent with the observation that, for small k, the lattice Uj\,(l) symmetry 
becomes a flavor non-diagonal axial symmetry in the continuum limit. 
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The hamiltonian strong-coupling analysis of Brower and Svetitsky ]Tj| was done 
for domain-wall fermions with a continuous fifth coordinate (05 = 0) in the limit 
L 5 — * 00, where L 5 is the size of the fifth dimension. We recall that in this double 



limit domain- wall fermions are expected to reduce to overlap-Dirac fermions [ ID ] . This 
work is not limited to a small hopping parameter. At zeroth order in the expansion 
in 1/g 2 they find a massive phase for M < 3 and a gap-less phase for M > 3. This 
leading-order result supports the mean- field argument of Sect. 6 that there is only 
one critical point at g — 00. (At that order there are no spatial couplings, and for 
each site the problem reduces to a free one-dimensional hamiltonian acting on the 



fifth coordinate. While the critical value in ref. |y| is M = 3, rather than M — 4, 
this is merely a technical difference stemming from the fact that time is taken to be 
continuous in ref. [fiT 



The order-l/g 2 result of ref. |TS] is summarized, for M > 3, by an effective 
low-energy hamiltonian in 3 + 1 dimensions which contains two terms, and (in their 
notation) reads Hf s + Hf te . The symmetry of H* a is that of naive fermions. Since, 
at that order, H^ s has only nearest-neighbor couplings, this can be explained by 
the fact that the O(k) term in the overlap-Dirac operator is proportional to the naive 
fermion one. At order 1/g 4 , H^_ s will contain also next-to- nearest neighbor couplings, 
and this may reduce the global symmetry to Ul(Ni) x Ur(N{). If we would ignore 
the second term, H^ e , this would be in agreement with our results. 

The second term in the effective hamiltonian, H^ e , explicitly breaks all axial sym- 
metries, having the same structure as an explicit mass term for the quarks. Although 
mathematically there is no direct conflict between this (hamiltonian) result and our 
(euclidean) result, physically the two results seem to be in conflict if both g and M 
are large (and therefore both results should be valid). We hope to resolve this issue 
in the future. To this end, it may be useful to carry out a euclidean strong-coupling 
analysis, for example using the method of ref. ||22j| . 

In conclusion, in this paper we have analyzed overlap-Dirac fermions with a small 
hopping parameter. While there are issues that require further work, a concrete 
picture of the phase diagram is beginning to emerge. In the future we hope to 
generalize the discussion to domain-wall fermions (with a non-zero and a finite 
L5), thus making closer contact with present-day numerical simulations. However, as 
we argued here, already in the limit where Liischer's chiral symmetry is exact, it is 
possible that a phase without massless quarks exist at fixed < M < 4 and large 
gauge coupling. 

Acknowledgements. We wish to thank Mike Creutz, Ikuo Ichinose, Ben Svetitsky and 
Marvin Weinstein for helpful discussions. MG would like to thank the Institute for 
Nuclear Theory at the University of Washington for its hospitality. This research is 
supported by the United-States - Israel Binational Science Foundation, and Y.S. is 
supported in part by the Israel Science Foundation, while MG is supported in part 
by the US Department of Energy. 

Appendix. Exponential localization. The Legendre expansion of (if 2 ) -1 / 2 (cf. eq. (|3])) 
was used in ref. [0] to prove the exponential localization of D. This expansion is very 
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informative, and may be used to obtain exponential-localization bounds on functions 
of D as well, in particular on D(2 — D)^ 1 . The Legendre expansion is given by [T3| 

(H 2 y 1/2 = C£t k <t> k (z) . (14) 

k>0 

Here <fik are Legendre polynomials, normalized such that max 4>k{z) = 1 for — 1 < 
z < 1. With slight adaptation the other ingredients are defined as follows. Let 
vq = max(H 2 ), uq = min(H 2 ), where extremization is done over the entire gauge- 
field configuration space. We take a (positive and) small enough hopping parameter 
k so that Uq > 0. Then 

z= u + v -2H^ 

Vq - Uq 

W') = !*±* (16) 

2 Vq ~ Uq 



with < t < 1 and 



At \ 1/2 

c=[-^—) • (17) 

Vq - Uq) 



Since \B\ < 4 and (for A = -1) 

H 2 = 1 + 2kB + 0(k 2 ) , (18) 

we have w = 1 + 8k + 0(k 2 ), u = 1 - 8k + 0(k 2 ), t = 4k + 0(k 2 ) and c = 1 + O(k). 
We may then write (for A = — 1) 

_ D = (c - 1) + c £ t fc (MZ) - b Y <f> h - X {Z)) , (19) 

k>l 

where b — n/t and Y = J^/il^C^ — B, cf. eq. (|3|). We comment that the Legendre 
expansion may be set up using any < u < u and v > v . The "best" values defined 
above, Uq and Vq, guarantee that tf k = 0(1). This is a natural relation because, as 
described below, t effectively plays the role of a hopping parameter. 
We now turn to the operator 

D(2- D)- 1 = Y / (D/2) n . (20) 

In order to obtain a bound on the kernel corresponding to eq. ( p0|) it is useful to 
regard t as an independent parameter. Doing so, we obtain the expansion 

D{2-D)- 1 = Y,t k V k , (21) 

fc>0 

where V k is defined by substituting eq. (|19"D into eq. (|20|) and collecting all terms 
involving an explicit factor of t k . At order t k , the operator encountered in the expan- 
sion of D (eq. fll9D ) allows for at most 2k hoppings |13[]. The same statement applies 
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to the new kernels: V k (x,y) = for \x — y\ > 2k where \x — y\ = \ x ^ — (the 
"taxi-driver distance"). Therefore 

D(2-D)-\x,y)= £ t k V k (x,y). (22) 

2k>\x-y\ 

What is still needed is a bound on V> k . We first observe that 

\<f> k (Z)-bY(j> k ^(Z)\ < 1 + 86. (23) 
In view of this bound we are lead to consider 

E = d + c(l + 86) ]T t k , (24) 

k>\ 

where d — |c— 1|, as well as 

E(2 - E)- 1 = £ (E/2T = E **** ■ ( 25 ) 

n>l k>0 

Again £ k is defined by collecting all terms involving t k using eq. (p4|). The double 
geometric series is easily summed, giving 

i e£ + (c + 86c — d)t , , 

£2-^ = — i — — , 26 

v ; 2 - d - (2 + c + 86c - d)t ' v ; 

and an explicit expression for follows by re-expanding eq. (^) in powers of t. 

Following the construction we see that, thanks to inequality (|23|), each term in 
T>k is bounded by a corresponding term in £ k . The latter is obtained if every factor 
of 4> k (Z) — bY (fik-i(Z) is replaced by 1 + 86 and every factor of c — 1 by d. Using 
eq. fl22|) we thus obtain a bound (for \x — y\ > 2) 

| J D(2- J D)- 1 (x,y)| < J2 **** = 

2k>\x-y\ 

d s t l*-wl/2] + t( c + 86c - d)st l^-yl/ 2 - 1 ] 



2 - d - (2 + c + 86c - d)t 



(27) 



where s = (2 + c + 86c — d)t/ (2 — d) and [I] is the smallest integer n such that n > I. 
Recalling the value of t — 4k + 0(k 2 ) (cf. eq. (|16|)), the desired exponential bound on 
the kernel of D{2 — D)~ l is thus given by eq. ( |27"1) provided we choose k < /to, where 
kq is the smallest value where the denominator in eq. fl2"T|) is zero. The denominator 
is strictly positive for < k < k,q since both t and d are O(k). 
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